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Abstract 

Monochromatic photons could be produced in the annihilation or decay of dark 
matter particles. At high energies, the search for such line features in the cosmic 
gamma-ray spectrum is essentially background free because plausible astrophysical 
processes are not expected to produce such a signal. The observation of a gamma-ray 
line would hence be a 'smoking-gun' signature for dark matter, making the search 
for such signals particularly attractive. Among the different dark matter models pre- 
dicting gamma-ray lines, the local supersymmetric extension of the standard model 
with small i?-parity violation and gravitino LSP is of particular interest because it 
provides a framework where primordial nucleosynthesis, gravitino dark matter and 
thermal leptogenesis are naturally consistent. Using the two-years Fermi LAT data, 
we present a dedicated search for gamma-ray lines coming from dark matter annihi- 
lation or decay in the Galactic halo. Taking into account the full detector response, 
and using a binned profile likelihood method, we search for significant line features 
in the energy spectrum of the diffuse fiux observed in different regions of the sky. No 
evidence for a line signal at the 5a level is found for photon energies between 1 and 
300 GeV, and conservative bounds on dark matter decay rates and annihilation cross 
sections are presented. Implications for gravitino dark matter in presence of small 
-R-parity violation are discussed, as well as the impact of our results on the prospect 
for seeing long-lived neutralinos or staus at the LHC. 



1 Introduction 



Gravitational evidences for dark matter exist from galactic to cosmological scales, but its 
particle nature still remains unknown [ll-Isj. The most popular dark matter candidate, a 
weakly interacting massive particle (WIMP), could in principle be detected directly through 
its scattering on atomic nuclei, indirectly by observing traces of its annihilation products 
in the cosmic-ray fluxes, or by its direct production at colliders. Since the stability of 
dark matter is by no means established, there is also the possibility to observe cosmic-ray 
contributions stemming from the decay of dark matter particles. However, despite lots of 
efforts, no unambiguous non-gravitational dark matter signal has been found so far. 

Indirect searches for dark matter annihilation or decay products in the cosmic gamma- 
ray, anti-particle and neutrino fluxes are limited not only by the instrument sensitivities, 
but also by our understanding of the astrophysical foregrounds. One of the few dark 
matter signatures that would unambiguously stand out of the astrophysical foreground 
are the monochromatic photons potentially produced in two-body annihilation or decay 
of dark matter. They would appear as a line feature in the otherwise continuous gamma- 
ray energy spectrum The observation of such a gamma-ray line would thus be of 
paramount interest for the understanding of dark matter in the Universe. 

The production of gamma-ray lines is expected in many theoretical dark matter models. 
The corresponding branching ratios are often very small, but this could be compensated by 
the high experimental sensitivity to such signals. In the popular scenario where dark matter 
is the lightest neutralino Xi of the minimal supersymmetric standard model (MSSM) [3], 
gamma-ray lines from y°y° 77 or XiXi ~^ 7-^° annihilations are only produced at 
the one-loop level |8-14|, and they would be mostly out of reach for current experiments. 



However, their production can be enhanced in non- minimal variations of the MSSM |l5l.ll6 
Alternatively, a large number of viable WIMP models exist, some of them predicting a 



sizable annihilation into monochromatic photons, e. g. si nglet dark matter 17|, hidden U(l) 



dark matter 1181-20 , effective dark matter scenarios 21 



into the higgs boson and a photon [22 



dark matter in models with universal extra dimensions 24, 2a 



scenarios with strong annihilation 
inert higgs doublet dark matter (23| . or Kaluza-Klein 

Additionally, prominent 



spectral features in the gamma-ray flux can be generated by the internal Bremsstrahlung 



or final state radiation potentially produced in dark matter annihilation |26H29 



An interesting example of decaying dark matter is the gravitino 1P3/2, which appears in 
locally supersymmetric extensions of the Standard Model. In scenarios where i?-parity and 
the lepton number are slightly violated and the gravitino is the lightest superparticle (LSP), 
thermal leptogenesis, gravitino dark matter and primordial nucleosynthesis are naturally 
Within this framework, the gravitino is no longer stable and decays on 
thus rendering the imprints of its decay potentially observable 
Interestingly, the line feature generated by the two-body 



consistent [30 1 



32 



cosmological time scales 
in the cosmic-ray fluxes 
decay into neutrinos and photons, 1P3/2 — ?■ 71/, is prominent for a wide range of gravitino 



-43 



masses 



33 



suggesting that the first observation of gravitino dark matter could happen 
via gamma-ray lines. 

In general, dark matter lifetimes of the order of 10^^-lO^^s, which is in the ballpark 
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of what is accessible by cosmic-ray experiments in the GeV-TeV regime, are obtained 
when the symmetry responsible for the dark matter stability is violated by dimension six 



operators generated close to the grand unification scale |44j . |45| . A typical example, which 



has a large branching fraction into monochromatic photons, is the hidden SU(2) vector dark 



matter [46j. Note that even if dark matter only decays into charged leptons at tree level 



like in the models designed to explain the PAMELA/Fermi LAT measurement s |47H 49 
gamma-ray lines produced at the one-loop level could potentially be observable [50| . 

A dedicated search for Galactic gamma-ray lines with energies between 30 and 200 GeV 
using the Fermi Large Area Telescope (LAT) data 51| has been presented in Ref. 52|, where 
the all-sky averaged diffuse gamma-ray flux measured between 7 Aug 2008 and 21 Jul 2009 
was considered. Searching for monochromatic features in the energy spectrum, their null 
result enable them to put constraints on the dark matter annihilation cross section and 
decay rate into gamma-ray lines. A similar analysis using the EGREtEL data has been 



performed for gamma-ray line energies between 0.1 and 10 GeV in Ref. [53 



The extra-galactic gamma-ray background (EGBG), as derived from the Fermi LAT 



data in Ref. 5J|, has also been used in Ref. 55|to put constraints on gamma-ray lines 



571 ]. In that case, the redshifted gamma- 



from cosmological dark matter annihilation [56|, 
ray lines lead to less pronounced features in the energy spectrum than the Galactic ones. 
The lack of spectral features in the EGBG was used to put limits on the corresponding 
annihilation cross section. Further gamma-line searches, coming from different observa- 
tions, and ranging from energies of 10 keV up to 10 TeV, were presented for annihilation in 
Ref. 58j and for decay in Ref. 



59| . General strategies for the detection of gamma-ray lines 



from dark matter annihilation with Fermi LAT and future experiments were discussed in 

Refs. leol-lei. 



The purpose of this paper two-fold: First, we will update and extend the gamma-ray 
line analysis of the Fermi LAT data presented in Ref. 52| to a larger energy range of 



1-300 GeV, searching for significant line signals that might come from dark matter anni- 
hilating or decaying in the Galactic dark matter halo, and derive conservative constraints. 
Although highly relevant for the gravitino dark matter scenario, no dedicated search for 
gamma-ray lines from decaying dark matter in the whole (9(GeV-TeV) energy regime has 
been performed so far. Furthermore, the large statistics of the Fermi LAT data allows us 
to further probe the low energy regime, i.e. energies from 1 to 10 GeV, and thus strength- 



ening the previous EGRET bounds on annihilating dark matter [53|. Such energies are 



of particular interest in the view of the recent hints from direct dark matter detection 



reported by the CDMS [63|] , CoGeNT [6J] and DAMA (65|, |66|] collaborations. Second, we 
will discuss the impact of our results on decaying gravitino dark matter. Covering previ- 
ously unexplored values of the gravitino mass, we will use our limits on gamma-ray lines 
to constrain the size of -R-parity violation, which can be translated into lower limits on the 
decay lengths of neutralino or stau NLSPs. 

The paper is organized as follows: In section 2, we present our analysis of the Fermi LAT 



^See http : //heasarc . gsf c . nasa . gov/ docs/ cgro/ cgro/ egret . htm"l[ 



3 



data. After briefly discussing the gamma-ray flux produced in the dark matter annihilation 
or decay, we present the strategy for data extraction and line searches in section 2.1. Section 
2.2 is devoted to a presentation of the resulting bounds on annihilation cross section and 
decay rate, whereas in section 2.3 we discuss details and caveats of our analysis. In section 3, 
we apply our results to gravitino dark matter and discuss prospects for NLSP observations 
at particle colliders. We then conclude in section 4. 



2 Gamma-ray line constraints from Fermi LAT 



Our line search is based on the measurements of the cosmic gamma-ray flux performed 
since August 2008 by the Large Area Telescope (LAT), the main instrument on the Fermi 
Gamma-ray Space Telescope. The LAT is a pair conversion detector, designed to observe 
gamma rays with energies ranging from 30 MeV to more than 300 GeV with an energy 
resolution of about 10%. It has a large field-of-view that spans around 20% of the sky, 
allowing an effective measurement of the all-sky diffuse fluxes. From the four event classes 
presented in the available public data, we only consider the recently released 'DataClean' 
event class, which provides the best rejection of charged cosmic rays. This is especially 
important when analyzing the diffuse gamma-ray flux at high energies. 

Gamma-ray lines may find their origin in the decay or annihilation of dark matter 
particles ^ in the Galactic halo. The observed gamma-ray flux at Earth is given by a line- 
of-sight integral depending linearly on the dark matter distribution Pdmi^) when considering 
dark matter decays 
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while in the annihilation case, it depends quadratically on pdm(^) 
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Here, is the dark matter mass, while T^^.^^ and (av)^^^^^ denote the dark matter decay 
width and annihilation cross section into gamma-ray lines, respectively. The coordinate s 
runs along the line of sight, in a direction spanning an angle ^ with respect to the Galactic 
center direction. The distance to the Galactic center r is related to the distance to the Sun 
s through r(s,^) = -\/ (rn — s cos + (s sin i^)^, where tq = 8.5 kpc is the distance of the 



Sun to the Galactic center 67|, l68| 



We consider the Navarro-Frenk- White (NFW) profile 52 
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Target region 


Geometry 


f. „ (in - 

J Ail P 

[GeV cm~^ sr] 


f. „ (in J^''""- 

J AS 2 P 

[GeV^ cm~^ sr] 


Halo 


|6| > 10° 


2.2 X 10^3 


8.3 X 10^2 


Center 


|/| < 36° and 5° < \b\ < 36° 
plus \l\ < 7° and \b\ < 5° 


6.9 X 10^2 


9.2 X 10^2 



Table 1: Geometry of our target regions in terms of Galactic longitude I and latitude b. 
The line-of-sight integrals Jp""^' and Jp""- integrated over the target regions are also shown, 
assuming a NFW profile. 



which is defined by {a, /3, 7) = (1, 3, 1) and = 20 kpc, as a reference for the dark matter 
distribution. For the sake of comparison, we also consider the isothermal profile with 
(a, /3, 7) = (2, 2, 0) and = 3.5 kpc, as well as the Einasto profile 70l-l72 



Pdm[r) = Poexp 







a 





with a = 0.17 



20 kpc 



(4) 



which is favored by the latest A^-body simulationsJl All profiles are normalized to 



To) 

Tann. 



Pdm[r 



0.4 GeV cm ^ at Sun's position [7J,|75j. The resulting line-of-sight integrals J^^^' and 



integrated over our target regions (see below), are summarized in Tab.[T]for the NFW 
profile. Note that substructures in the Galactic dark matter distribution can change the 
angular profile and magnitude of the gamma-ray emission (see e.g. Refs. 62, 72||). If this 
is the case, our results can straightforwardly be generalized by adopting the corresponding 



1 n ydcc/ann. 

values for Jp 



2.1 Methods 

The gamma-ray events that enter our analysis are selected from the 'DataClean' event 
class of the Fermi LAT data measured between 4 Aug 2008 and 17 Nov 201o|fl From all 
events recorded by the Fermi LAT, we select energies between 100 MeV and 500 GeV, and 
apply the zenith angle criterion 6 < 105° in order to avoid contamination by the Earth's 
Albedofl Below 1 GeV, the fast variation of the effective area as well as the presence of 
small and sharp jumps in the energy spectrum (presumably related to side effects of the 



energy reconstruction [76|, |77|) introduce spurious effects in our line search. Furthermore, 
the nominal energy range of the Fermi LAT ends at 300 GeV jsij. As a consequence, we 
concentrate on gamma-ray line energies between 1 and 300 GeV. Systematic uncertainties, 



^However, see Ref. [73[ for observational arguments in favor of cored dark matter profiles. 

•^The event data, as well as the corresponding information about the instrument response functions 



concerning angular and energy resolution, are available at (iiittp : //f ermi . gsf c . nasa . gov/ ssc/ data/ 



^ These selections are made using the Fermi Science Tools v9rl8p6. 



see 



http://feriiii.gsfc.nasa.gov/ssc/data/analysis/software/ As cut in gtmktime we took 



DATA_QUAL==1 && LAT_C0NFIG==1 && ABS(R0CK_ANGLE)<52. 
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10° 10^ lO' 10° 10^ lO' 

E [GeV] _B., [GeV] 

Figure 1: Left panel: Observed counts in the halo and the center regions, binned into 50 bins 
per energy decade (200 bins per decade are used in our hkelihood analysis). Right panel: 
Fermi LAT energy resolution AE/Ej, in terms of FWHM and 68% containment (t^e^ 
derived from the IRF of the 'DataClean ' event class. 



as well as limitations of the analysis method, are discussed in section 12.31 in light of the 
data. 

When looking for dark matter signals in gamma-ray maps that cover the whole sky, it 
is critical to choose a target region which maximizes the corresponding signal-to-noise ratio 
S/N. The two regions that we found giving a very good S/N for decaying dark matter 
('halo region') or annihilating dark matter ('center region') are summarized in Tab. [H 
Although optimized for the NFW profile, they also yield good S/N for the Einasto and 
isothermal profiles. 

Since the selection procedure leaves us with a large number of 1.3 x 10^ (5 x 10^) events 
above 1 GeV for the halo (center) region, we perform a binned analysis of the data. To 
this end, we distribute the events into 200 logarithmically equally spaced energy bins per 
decade, and sum over the angles. This gives a sequence of count numbers q G Nq, which is 
illustrated in the left panel of Fig. [1] for both target regions. Note that we do not perform 
a point source subtraction in this work. A proper treatment would mask out only (9(5%) 



of the events [52[, and hence only marginally affect our results. 

The spectral feature produced by a gamma-ray line can be inferred from the Fermi 
LAT instrument response function (IRF). Its most recent version, 'Pass6 version 3', was 
determined using Monte Carlo generated samples of photon events between 18 MeV and 



562 GeV, and includes effects measured in-flight, see Refs. [78|, |79|. It contains the point- 
spread-function (PSF), as well as the energy dispersion V{E,E^) which describes the dis- 
tribution of the reconstructed energies E as a function of the physical photon energy E.y. 
In order to integrate out the implicit dependence of the energy dispersion on the event im- 
pact angle with respect to the detector axis, V^E, E^) is averaged over this impact angle 
weighted by its distribution in our data sample. The resulting full width at half maximum 
(FWHM) of the energy dispersion is shown in the right panel of Fig. [H where the 68% 
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E [GeV] E [GeV] 

Figure 2: Left panel: Integrated exposures £int for both target regions, following the weighting 
method described in the text. Right panel: Differential gamma-ray flux as observed in the 
two target regions. 



containment energy resolution (7^^° ^^^^ shownj^l Lastly, since our two samples have a 
large sky coverage, the effects of the PSF, which features a 68%-containment angle well 
below 1° at energies above 1 GeV, can be safely neglected. 

The profile likelihood method 80|, |81[ is used to calculate the significance of a potential 
gamma-ray line contribution to the observed gamma-ray flux. The likelihood function is 
defined as 



L{n\c) = Ui=i^^,„^i^Pf,Xci) , (5) 

where, for each energy bin i, /ij and q denote the expected and observed count numbers, 
respectively, and P^{k) = fi^e~^/k\ is the Poisson probability mass function with mean 
value fi. The data are modeled by a simple power law with a spectral index —7 and 
a normalization /3, plus a line signal at fixed energy i?^ with the normalization a. The 
power law background is supposed to take into account not only the gamma-ray flux of 
astrophysical origin, but also the gamma-ray continuum which might be produced in the 
dark matter annihilation or decay additionally to the sharp gamma-ray line. Since the 
power law is only locally a good approximation to the background fluxes, we use a sliding 
energy window in the fitting procedure. For gamma-ray energies < 10 GeV we only 
take into account the energy bins i G [io?"^!] th.ai lie within ±2cr^^^° around E^, whereas, 
due to limited statistics, at high energies E^ > 100 GeV bins within a wider energy window 
from roughly ^E^ to 3E^ are considered. In the range 10 GeV ^ E^ < 100 GeV the energy 
window size is interpolated between these two cases. Furthermore, events above 300 GeV 
are always excluded from the analysis, which leads to energy windows asymmetric around 

^Our resulting FWHM is larger than the one discussed in Ref. [H^l by a factor of 1.1-1.4, which 
presumably affects our limits at the level of 5%-20%. The difference is probably due to different energy 
reconstruction methods underlying the public 'DataClean' event class used in this work and the dedicated 
analysis in Ref. jslj. However, the 68% containment agrees well with Ref. [Hj. 
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the highest gamma-ray hne energies. The expected number of counts in each energy bin 
is given by 

/i,(a,/3,7)= / ' dE' {P£UE')E'"^ + a£UE,)V{E\E,)) , (6) 
Je- 

where the Ef denote the boundaries of each energy bin z, £^int is the energy dependent 
integrated total exposure, and we neglect the effects of the energy dispersion on the power- 
law background flux. To obtain ^int, we integrate the exposure maps, which we derived 
with the Fermi Science Tools@ over the observed region of the sky, weighted by the angular 
profile of the expected gamma-ray line signal from dark matter decaying or annihilating 
inside the Galactic halo. The resulting integrated exposures are practically identical for 
decaying and annihilating dark matter signals, and are shown in the left panel of Fig. [2] for 
the halo and center regions. They can be approximated by £\t^t{E) — {£{E))AQ, where 
the target region averaged exposures are in the range {S{E)) ~ 5-8 x 10^° cm^ s, and 
AQ = 10.4 sr (1.30 sr) is the solid angle spanned by the halo (center) region. The observed 
differential gamma-ray fluxes for both target regions, as inferred from the exposure and 
count maps, are shown in the right panel of Fig. [2j 

Depending on whether the normalization a of the gamma-ray line is kept free or fixed, 
the likelihood function in Eq. must be maximized with respect to three or two param- 
eters, respectively: 

L(c) = sup{L(At|c) : a, /3, 7 e M+} or L„(c) = sup{L(m|c) : /3, 7 G M+} . (7) 

Defining Aq,(c) = La{c)/L{c), the 95% CL. interval upper limit of the gamma-ray line 
normalization is given by the largest value of a which satisfies — 21ogAQ(c) < 4. The 
parameter a is defined such that it can be directly identified with the corresponding gamma- 
ray line flux, i.e. = a. In the absence of a gamma-ray line signal at the 5(j level, 
which would correspond to — 21ogAa=o(c) > 23.7, we quote the upper limit of the 95% 
CL. interval as the upper limit on the gamma-ray line flux. Note that the — 21ogAQ=o 
value corresponding to a 5(j-confidence level slightly deviates from the naively expected 25, 
because the normalization of the gamma-ray line is bounded from below, i.e. a> 0. The 
adopted value can be derived by recognizing that in absence of a gamma-ray line in 50% 
of the cases the optimal fit would actually require an unphysical negative normalization 
of the gamma-ray line, since for large enough count numbers statistical fluctuations are 
symmetric around their mean. These cases yield — 21ogAQ,=o = 0, whereas in the other 
50% of the cases, the — 21ogAQ=o would follow a conventional x^-distribution with one 
degree of freedomQ 



^In gtltcube we use a zenith cut zmax=105 to properly calculate the exposure also in cases where the 
field-of-view of the satellite intersects with the Earth's Albedo. However, due to our above cut on the 
satellite rocking angle (see footnote H]) this has only a small effect on the total exposure. 

^We confirmed this behavior with a Monte Carlo analysis. Details about the likelihood-ratio and 
confidence intervals can be found e.g. in Refs. [sH, [s^]. Note that wc do not use a "trial factor" as 



discussed in Refs. |76l. l84l| . which only becomes relevant when constraining models that predict a multitude 



of lines or when claiming a line discovery. 
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Upper limits on gamma-ray line flux at 95% and 99.7% C.L. 



Halo region 
Center region 




[GeV] 



Figure 3: Our results for the 95% and 99.7% C.L. upper limits on the monochromatic photon 
Rux from the halo and the center regions, as function of the gamma-ray line energy Ej . Note 
that due to instrumental effects these limits are expected to be too strong at lower energies 
in some cases (cf section [0]) . The corresponding conservative limits are show in Tab. [3| (as 
well as in the Figs. [Hand I^)- 



The likelihood analysis is performed using PyFITSlfl and PyMinuit, an interface to the 
function minimizer MINUIT 1.7.9j^ 



2.2 Results 

No gamma-ray lines with 5a significance are found in the 1-300 GeV energy range. Details 
and caveats of the analysis are discussed in section [2731 Our results for the 95% and 99.7% 
C.L. limits on the gamma-ray line flux are shown in Fig. |3]as a function of the gamma-ray 
line energy for both halo and center regions. As discussed in section 12.31 instrumental 
effects forbid us to take these limits at face value; the conservative limits summarized in 
Tab. [3] should be used instead. These are derived by quoting the weakest limit we obtain in 
given energy bands. Together with Eqs. ([T]) and (l2l), it is then straightforward to derive the 
corresponding constraints on the dark matter decay width or annihilation cross section. 

Limits on the inverse dark matter decay width into monochromatic photons and neutri- 
nos, T^\^^^, as derived from the halo region observations, are presented in Fig. HI assuming 
a NFW dark matter profile. The gray line shows the bounds as they directly follow from the 
flux limits in Fig. [3l In most of the considered energy range, they fluctuate around 10^^ s. 



*PyFITS is a product of the Space Telescope Science Institute, which is operated by AURA for NASA, 
see http : //www. stsci . edu/resources/ sof tware_hardware/pyf its 



See|http : / / code . google . com/p/ pyminuit/ and http : / / seal . web . cern . ch/ seal/ snapshot/work-packages/mathld 



Target region 


pprofile ;p 




ijjtjj^'y'y 




Iso. NFW Bin. 


Iso. NFW 


Ein. 


Halo 
Center 


0.96 1.00 1.00 
1.03 1.00 0.91 


1.03 1.00 
1.65 1.00 


0.89 
0.64 



Table 2: Rescaling factors for tlie bounds on the dark matter decay width and annihilation 
cross section (shown in Tab. 0), corresponding to the different halo profiles. In most cases, 
the weakest bounds are obtained for the isothermal profile (Iso.), the strongest bounds for 
the Einasto profile (Ein.), whereas the NFW profile leads to intermediate constraints. 



The slight improvement with respect to the previous Fermi LAT hmits [52] (dashed line) is 
due to the increased statistics after two years of data taking. The black dots represent the 
weakest hmits obtained when varying the gamma-ray hne energy within different adopted 
energy bands. Except at the highest energies, these values fluctuate around 6 x 10^*^ s, and 
are summarized in Tab. [31 Adopting halo profiles different from the NFW profile would 
change the limits by a few percent. The corresponding rescaling factors for the Einasto 
and the isothermal profiles are summarized in Tab. |2l and their impact on the limits is 
illustrated by the blue band in Fig. HI Note that the normalization of the dark matter 
profile at Sun's position is another source of uncertainty, which is not captured in the 
rescaling parameters of Tab. El We do not plot the somewhat weaker limits on the decay 
width obtained from the center region, which can be found in Tab. |3lfor completeness. 

In Fig. |5lwe present the limits on the dark matter annihilation cross section into gamma- 
ray pairs, {av)^^^^^, as derived from the center region observations. As ihustrated by 
the blue band, the uncertainties coming from the dark matter density profile are much 
larger than in the case of dark matter decay and of the order of ~ 50%, cf. Tab. [21 For 



comparison, we also show the previous limits obtained by EGRET 53|o The present 
analysis improves these limits by up to an order of magnitude, which is essentially due 
to the increased statistics of Fermi LAT with respect to EGRET. Furthermore, at high 



energies, we slightly strengthen the previous Fermi LAT constraints [5^. Our results are 
summarized in Tab. [31 It is interesting to note that, assuming the NFW profile, the limits 
on the dark matter annihilation cross section derived from the halo region are only a factor 
~ 2 weaker than the limits derived from the center region. 

Using kinematic considerations, our results can be straightforwardly translated into 
constraints on two-body decay or annihilation into a monochromatic photon plus a mas- 
sive neutral particle N with mass rriisf. This is relevant for the common case of annihila- 
tion/decay into 7^". However, in principle rriN is a free parameter, and if rriN is close to the 
dark matter mass , the line signal can be even strongly enhanced in some scenarios 5^ ■ 



-'^^Note that the cross section hmits shown in Fig. 8 of Ref. [53[ are incorrectly too strong by a factor of 
ten, whereas the flux limits presented in Fig. 7 are correct [priv. comm. with the authors]. The EGRET 
cross section limits that we present in Fig. [5] are derived from Fig. 7 of Ref. [ssj . 
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Lower limits on DM inverse decay width into (halo region) 




[GeV] 

Figure 4: Lower bounds on the dark matter inverse decay width into monoclu-omatic photons 
and neutrinos, T^^^^, as a function of the dark matter mass m^, derived from the halo region 
fluxes assuming the NFW dark matter profile. The gray-solid line shows the 95% C.L. limits 
as directly derived from the line flux limits shown in Fig. O The black dots show the weakest 
limits obtained in the adopted energy bands and are listed in Tab. O For comparison, the 
previous Fermi LAT limits from Ref. jsl] are also shown by the red-dashed line (rcscalcd to 
our decay channel). The barely visible blue bands illustrate how the bounds change when 
using the isothermal or Einasto dark matter profiles instead. 



Defining the quantity 



x = i + di + ^, (8) 

which is bounded from below by x > 2, our limits for the decay into translate into 
limits on the decay into via 

X 

r^^^Ar(m^) = -r^^^^{E^) , where = x . (9) 
For the annihilation scenario the cross section limits on ipip — )■ are obtained by 

X 

{av)^^_,^N{m^) = y {crv)^^^^^{E^) , where m^ = -E^ . (10) 



Models that are partially constrained by our limits include singlet dark matter 17| and 

Scenarios that remain practically 



SU(2) vector dark matter [46|, but also Refs. [I9|-l22, 
unconstrained include the standard MSSM scenario 
particular case of gravitino dark matter will be discussed in the next section. 



but also Refs. 15, 16, 23-25 . The 
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Upper limits on DM annihilation cross section into 77 (center region) 
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Figure 5: Upper bounds on the annihilation cross section into gamma-pairs, {av)^^^^^, as 
a function of the dark matter mass m^, derived from the center region fluxes assuming the 
NFW dark matter profile. The gray-solid line shows the 95% CL. limits as directly derived 
from the line Rux limits shown in Fig. [3l The black dots show the weakest limits obtained in 
the adopted energy bands and are listed in Tab. O For comparison, the previous Fermi LAT 
limits from Ref [52] as well as the limits derived from EGRET observations of the Galactic 
center jS^] are also shown by the red-dashed and the black-dotted lines, respectively. The blue 
bands illustrate how the bounds change when using the isothermal or Einasto dark matter 
profiles instead. 



2.3 Discussion 

One crucial assumption underlying our analysis is that the background flux in the different 
considered energy windows can be well approximated by a power-law. This assumption 
is most likely to break down in cases where the statistics is very good. In order to check 
the validity of a power-law ansatz, we show in Fig. [6] the x^/d-o.f. of the background-only 
(green lines) and of the background-plus-signal (red lines) fits, as function of the gamma- 
ray line energy]^ The grey band corresponds to a p- value of > 5%. For the center region 
the fits are essentially in agreement with the data over the whole energy range. However, 
p-values significantly smaller than 5% occur at energies between 1 and 10 GeV (as well as 
at high energies close to 300 GeV) when considering the halo region, which has a three 
times larger statistics than the center region. Assuming that the astrophysical gamma- 
ray fluxes follow smooth bended power- laws, this tension points to an instrumental effect, 
presumably related to the energy reconstruction of gamma-ray events. 

^^The smallness of the differences between the x^/d.o.f. of the background-plus-signal and background- 
only fit at high energies comes from the fact that the values are actually dominated by the background 
and not by the narrow signal. 
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Center region 



/d.o.f. (BG+line fit) 
/d.o.f. (BG only fit) 



Significance \/-2iogA^ 




Halo region 



X^ /d.o.f. (BG+line fit) 
X^ /d.o.f. (BG only fit) 



Significance v'-21ogA„^o 




Figure 6: Significance of a potential line signal (black line) in the center (left panel) and halo 
(right panel) region in terms of standard deviations. We also show the /d.o.f. values of fits 
with (red line) and without (green line) line signal. The gray shaded area corresponds to a 
p-values of > 5%. Especially at low energies in the halo region, significant deviations occur 
from what is statistically expected in absence of a signal, as discussed in section [ 



The black lines in Fig. |6] show the significance (in terms of standard deviations) of 
line signals with respect to a power-law background. Over the whole energy range we are 
performing (9(50) statistically independent tests for a line. Hence, we expect to see a few 
cases with 2a significance, and practically no case with 3a significance. This expectation 
is marginally compatible with the observations of the center region, and is clearly violated 
for the halo region, where four > 3a excesses are observed below 20 GeV. However, an 
interpretation of these excesses in terms of gamma-ray lines would be suspicious: They 
occur at low energies in the halo region, where the statistics is best and hence the im- 
pact of instrumental effects most relevant, as indeed indicated by the bad x^/d.o.f. of the 
background-fits at this energy range. Moreover, most line-like features appear, although 
with a different significance, at the same position in both, the halo and the center region, 
indicating again a systematic effect. Without further understanding and reduction of in- 
strumental effects, no conclusive statement about the existence of gamma-ray lines (below 
~ 5(7 significance) in this energy range can be made. 

The limits presented in Figs. HI |5] and Tab. |3] are conservative in the sense that we 
always quote the weakest limit we obtain in given energy bands. Especially in the halo 
region at energies below 20 GeV, this ensures that our limits are dominated (and weakened) 
by the presence of the line-like artefacts in the energy spectrum. A future removal of these 
artefacts from the data is expected to lead in general to a strengthening of the presented 
limits. Note that the presence of a real gamma-ray line signal significantly stronger than 
the artefacts would have lead to limits dominated by the line signal itself. 

For two different gamma-ray line energies and for both center and halo regions, we 
show in Fig. [7] the line signals that are excluded at 95% C.L. (for which — 21og Aq(c) = 4, 
see above), and compare them with the best-fit background-only fluxes. The gamma-ray 
line energies are chosen such that the corresponding limits dominate in their respective 
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Figure 7: Line signals that are excluded at 95% C.L. (solid line), compared to best-Gt 
background-only model (dotted line) and data. We show results for the center (top pan- 
els) and halo (bottom panels) region. The gamma-ray line energies are = 1.5 GeV (left 
panels) and = 180 GeV (right panels). In the left panels, the signal-to-background ratio 
is 2.6% (top) and 1.5% (bottom). 



energy bands in case of the halo region. Furthermore, the energy ranges shown in the plots 
correspond to the energy windows used in the fitting procedure, and we rebinned the data 
for better visibility. At low energies the gamma-ray line constraints essentially come from 
the non-observation of a strong bending in the measured energy spectrum (with signal-to- 
background ratios as small as 1.5%), whereas at high energies, where the considered energy 
window is larger, the constraints come from the non-observation of a well defined bump. 

In summary, due to the apparent presence of line-like artefacts at the few percent level 
in the energy spectrum of the public 'DataClean' event class of the Fermi LAT, the flux 
limits from Fig. [3] are expected to be too strong in some cases, especially at energies below 
~ 20 GeV. Instead, the limits presented in Fig. HJ [S] and Tab. [3] should be used. 



Further more general uncertainties and systematic effects enter our analysis [76|, I77j : 
The systematic uncertainties in the effective area are energy dependent and increase from 
5% at 560 MeV to 20% at 10-100 GeV in case of the 'DataClean' event class Q (note that 
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[GeV] [ph/(cm2 s sr)] ^^\^u N (av)^,^,^^^ [cm^s ^ 

halo center halo center halo center 

xlO-9 xlO-9 xlO^s xlO^s X 10-28 X 10-28 



1.0-1.4 


11 


40 


6.4 


4.4 


0.25 


0.10 


1.4-1.9 


11 


60 


5.0 


2.3 


0.41 


0.24 


1.9-2.5 


7.9 


29 


4.7 


3.3 


0.63 


0.26 


2.5-3.5 


5.0 


18 


5.9 


4.6 


0.63 


0.22 


3.5-4.7 


3.9 


17 


5.9 


3.4 


0.82 


0.42 


4.7-6.5 


2.1 


8.8 


6.6 


4.7 


1.3 


0.49 


6.5-8.8 


2.4 


8.8 


5.2 


3.6 


1.7 


0.69 
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0.60 


3.1 


12 


6.6 


1.4 


0.60 


12-16 


0.98 


3.0 


5.2 


4.4 


4.2 


1.4 


16-22 


0.99 


4.4 


5.1 


2.9 


4.4 


2.2 


22-31 


0.42 


1.5 


6.6 


5.9 


6.2 


1.5 


31-42 


0.42 


1.7 


6.5 


3.0 


6.3 


5.1 


42-57 


0.21 


1.5 


7.1 


3.3 


11 


4.8 


57-77 


0.20 


0.92 


6.2 


3.0 


17 


9.4 


77-106 


0.18 


0.90 


5.0 


2.6 


29 


15 


106-144 


0.17 


0.78 


4.7 


2.6 


31 


18 


144-197 


0.14 


0.29 


3.3 


3.9 


72 


18 


197-300 


0.10 


0.59 


4.1 


1.2 


96 


91 



Table 3: Conservative upper limits on the Qux in gamma-ray lines J-y, as well as on the decay 
width r^_^^y and the annihilation cross section {av)^^^^-y of dark matter into monochro- 
matic photons, for line energies = 1-300 GeV. In each of the indicated energy bands, we 
varied to find the weakest possible 95% C.L. limit. The flux limits for all values of E^ are 
shown in Fig. O (see also section 2.3). The limits are evaluated for the halo and for the center 
region separately, and hold for the NEW profile. Other dark matter profiles just require a 
rescaling of the bounds as shown in Tab. 



these uncertainties amount mainly to an overall shift in the effective area). Uncertainties on 
the level of around 20% are expected to continue up to energies of hundreds of GeV |76l. |77| . 



but one has to keep in mind that a precise determination of the effective area at these high 
energies is still lacking due to limited statistics. However, even values reasonably larger 
than 20% would only mildly affect our conclusions. Over the whole considered energy 
range, the systematics of the LAT effective area consequently translate into uncertainties 
for our dark matter limits that range between 5% and ~ 20%. The absolute energy scale of 



the Fermi LAT is determined up to [48[. This translates into an additional uncertainty 



on the dark matter annihilation cross section {av)^^^^^ of ^20%' decay width 

T^^jiy of ^^0%, respectively. 

Finally, we have also checked that when extending our analysis to 500 GeV no significant 
line signals at high energies are revealed, and the limits below the highest energy band 
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presented in Tab. |3] are only marginally affected. Considering the 197-300 GeV (300- 
500 GeV) energy band, we obtain in this case > 6.0 x 10^^ s (1.7 x 10^*^ s) for the 

halo region, and (av)^^^^^ < 4.1 x 10"^'' cm'^ (1.54 x 10"^^ cm^ s"^) for the center 
region. However, we emphasize that energies above 300 GeV lie outside of the nominal 
energy range of the Fermi LAT. 



3 Consequences for Gravitino Dark Matter 



Local extensions of the supersymmetric standard model predict the existence of the gravi- 
tino, the supersymmetric partner of the graviton. Depending on the mechanism of super- 
symmetry breaking, it can be the lightest superparticle (LSP). In this case, it provides a 
natural dark matter candidate 85[ , and is an interesting alternative to the standard WIMP 
scenario. 

The dominant contribution to the gravitino production in the early Universe comes 
from the 2-to-2 QCD scatterings, which are taking place in the thermal bath at reheating. 
The resulting energy density is given by 



^3/2 



^ 100 GeV \ 

m3/2 J 



ITeV 



1010 GeV 



(11) 



where Tr is the reheating temperature, while 777.3/2 and are the gravitino and gluino 



masses, respectively. To leading order in the gauge couplings, C ~ 0.5 |86H88| j I Re- 



markably, the observed dark matter relic density VLduJi^ = 0.11 89| is obtained for typical 



supersymmetric parameters, 7723/2 ~ 100 GeV and r?7^ ~ 1 TeV, and for reheating tempera- 
tures Tr ~ 10^° TeV as required for thermal leptogenesis to generate the observed baryon 
asymmetry 
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However, the above scenario is not free of problems: If i?-parity is conserved, the decay 
of the next-to-lightest supersymmetric particle (NLSP) into gravitino and standard model 
particles is suppressed by the Planck scale. As a consequence, the NLSP becomes very 
long-hved and can dramatically affect the successful prediction of the standard big bang 
nucleosynthesis (BBN) scenario 92|-l98l]. Several alternative scenarios have been proposed 
to consistently accommodate the rma l leptogenesis, gravitino dark matter and BBN, like 
assuming sneutrino [9^ or stop 100 1 NLSPs, diluting the NLSP abundance by ent ropy 



production |lOl| . or introducing new NLSP decay channels into hidden sector states [102| . 
Il03i] . Also, the reheating temperature can be considerably reduced by considering non- 
thermal leptogenesis scenarios, as for example recently proposed in mechanisms where the 
universe is reheated th rough the decay of non-relativistic right-handed neutrinos produced 



in false vacuum decay |104 



105 



Alternatively, it has been proposed that /2-parity could be mildly violated (30|. Such 



a violation could be e.g. produced by dynamical breaking of 1]{1)b-l at low scales in the 

"'^^Note that C has 0(1) uncertainty due to unknown higher order contributions and nonperturbative 
effects ii. 
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hidden sector 106l |. It would lead to a rapid decay of the NLSP before the onset of BBN, 
thus rendering its impact on the standard cosmological picture negligible. Notably, even 
though the gravitino is not stable anymore in such a case, it still constitutes a viable dark 
matter candidate because its decay into standard model particles is doubly suppressed 
by the small -R-parity breaking parameter as well as by the Planck mass 3l|, leading to 
lifetimes longer than the age of the Universe. 

This i?-parity violating scenario is particularly attractive because it opens up the possi- 
bility to probe gravitino dark matter by searching for its decay products in the cosmic-ray 



fluxes. In addition to cosmic anti-matter (35|, l36| and neutrino (38j fluxes, the produced 



gamma-ray flux |30H34|. |36[ is of high interest because, besides the continuous component 
generated by the fragmentation of the Higgs and gauge bosons, the gamma-ray energy 
spectrum typically features an intense line, arising from the 'Ji^ two-body decay 

already at tree level. 

In the rest of this section, we will apply the above Fermi LAT limits on gamma-ray 
lines to the decaying gravitino scenario. First, we will summarize the relevant aspects 
of bilinear i?-parity violation, which will form the theoretical framework of our analysis. 
Second, we will present limits on the size of i?-parity violation, and discuss prospect for 
see ing long-lived neutralino and stau NLSPs at the LHC, extending the analysis made 
in ITotI. 



3.1 i?-parity breaking model 

The supersymmetric standard model with bilinear i?-parity breaking is, in addition to 
the standard SU(3)cXSU(2)x,xU(l)y gauge interactions, defined by mass-mixing terms 
between lepton and Higgs fields in the superpotential 

W = PFmssm + fiiHJi , (12) 

as well as by the mass-mixing terms in the scalar potential as induced by supersymmetry 
breaking 

C = C^f^ + BM + m^, l\Hd + h.c. . (13) 

Here, VFmssm and -^^^1^^ are the usual -R-parity conserving MSSM superpotential and 
scalar Lagrangian, Hu/d are the up/down-type Higgs doublets, U the lept on d oublets, and 



lii, Bi and mf^ are the i?-parity violating couplings. As proposed in Ref. [lQ7|, it is conve- 
nient to work in a lepton-higgs basis where the mass mixings /ij, Bi and m^^ are traded for 
i?-parity breaking Yukawa couplings. This can be achieved by a supersymmetric rotation 
of the superfields, followed by a non-supersymmetric rotation of the scalar fields alone. 
Besides the i?-parity breaking Yukawa couplings, which will become relevant for the decay 
of the stau NLSP below, these rotations also generate mass-mixing terms between the neu- 
tralinos and the neutrinos. This neutrino-neutralino mixing, which can be parameterized 
by the dimensionless parameter finally induces the gravitino decay mode 11)^/2 7^ as 



17 



well as the decay of neutralino NLSPs (see Ref. [io3 for a definition of C in terms of the 
bilinear i?-parity violating couplings /ij, Bi and m^^). 

We will consider two typical sets of boundary conditions for the supersymmetry break- 
ing parameters of the MSSM at the grand unification (GUT) scale for definiteness. The 
first one corresponds to equal scalar and gaugino masses 

(A) mo = mi/2, flo = 0, tan/3 = 10 , (14) 

for which the bino-like neutralino Xi is the NLSP. In the second one, which corresponds to 
no-scale models or gaugino mediation, 

(B) mo = 0, mi/2, cto = 0, tan/3 = 10, (15) 

the lightest stau fi is the NLSP. In both cases, tan /3 = 10 has been chosen as a repre- 
sentative value, and the trilinear scalar coupling has been set to zero for simplicity. 
The universal gaugino mass mi/2 remains as the only independent variable. For both sets 
of boundary conditions, the gaugino masses Mi_2,3 satisfy the following relations at the 
electroweak scale 

Mo 

^^5.9, 7^=^1.9. (16) 



Electroweak precision tests (EWPT) yield important lower bounds on the superparticle 



mass spectrum [88|]. For a neutralino NLSP, the most stringent constraint comes from the 



Higgs potential. The universal gaugino mass mi/2 is required to be high enough in order 

This implies the 



108 



for the Higgs mass to fulfills the LEP lower bound m^ > 114.4 GeV ^ ^ 

lower limit m^o > ISOGeVS However, allowing negative ao or scalar masses much larger 
than mi/2 at the GUT scale would weaken this limit, and we will take m^o > 100 GeV 
as a lower bound for the neutralino mass subsequently. In the stau NLSP case, the lower 
bound comes from th e abs ence of pair production of heavy charged particles at LEP and 
reads ruf, > 100 GeV fioi ]. 

For a typical effective neutrino mass mi = 10~^eV, successful thermal leptogenesis 

lO^GeV |9l|. Using Eq. together 



109| , this implies a lower bound for 



requires a minimal reheating temperature of Tr 
with a lower bound on the gluino mass rrig > 500 GeV 
the gravitino mass m3/2 > 5 GeV. 

In addition to the lower limits, NLSP mass upper limits follow from the requirement 
that the gravitino does not overdose the Universe. Indeed by rewriting Eq. ( ITT]) one obtains 
the constraint 



mNLSP — 310 GeV 



m3/2 / lO^ GeV y/' 

0.2 J V 100 GeV j [ Tr J 



(17) 
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Note that m^o ~ Mi with good accuracy [107 
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where ^ = mNLSp/''^g is implicitly fixed by the supersymmetry breaking boundary condi- 



1/2 

tions [88[. Since jtinlsp requiring the gravitino to be the LSP induces absolute 



''3/2' 

upper bounds on the NLSP masses. In the case of the neutralino NLSP, Eq. ( ITTl) implies 
^x? ~ 690 GeV for ^ = 1/5.9, and is essentially independent of mo and tan/3. For the 
stau NLSP, tan/3 = 10 yields C, = 1/6.2, which consequently leads to the more stringent 



bounds ruf^ < 615 GeV. Note that there is a strong dependence on tan/3 in that case 88l |. 
and that decreases with increasing tan /3. 

3.2 Fermi LAT constraints on the NLSP decay length 

The gravitino decay mode relevant for the discussion in this paper is the decay into pho- 
tons and neutrinos via -1/^3/2 — )■ 7Z/. However, decays into W^l^, Z^p and h?p final states 
are also possible, and have often large branching ratios. Below the corresponding kine- 
matic thresholds, three-b ody decay with intermediate massive gauge bosons can become 



important in some cases (42 



Additionally to the gamma-ray line from 11)^12 — )■ 7z/, the different decay modes will 
produce a gamma-ray continuum, coming from Higgs and gauge boson fragmentation as 
well as from final state radiation of the charged leptons. This continuous spectrum has a 
cutoff at the position of the gamma-ray line, namely at half of the gravitino mass. Since we 
required in the above line search that all gamma-ray fluxes except the line follow a power- 
law locally around the line's energy (c/. section [2?T]) . a too large continuum contribution 
from gravitino decay would render our line search strategy invalid. We expect this to be- 
come relevant for branching ratios into gamma-ray lines below 10""^ . . . 10^^, and postpone 
a detailed analysis of this effect to future work. Using the branching ratios presented in 



Refs. [36|, 138|, we find that this problem can become severe for gravitino masses above a 
few hundred GeV, ^3/2 ^ 300 GeV. However, one has to keep in mind that the exact 
branching ratio into lines is in principle mo del- dependent. In what follows we assume for 
simplicity that no continuum contribution disturbs our line search, and we present limits 
over the full range of accessible gravitino masses. 

The gravitino inverse decay rate into photon/neutrino pairs is given b}{3 31, 107 



32V2GfM| MfMf 



where a is the electromagnetic fine structure constant, Mp = 2.4 x 10^^ GeV the reduced 
Planck mass, and Gp = 1.16 x 10~^ GeV~^ is the Fermi constant. Using this expression, 
the limits on the gamma-ray line flux can be translated into upper limits on the i?-parity 
breaking parameter (. A conservative bound on ( is obtained for a given gravitino mass 
when using the corresponding maximally allowed bino mass. The latter is obtained by 
combining Eq. (ITT]) and ( IT6|) considering the lowest reheating temperatures allowed in the 
thermal leptogenesis scenario, i.e. Tji ~ 10^ GeV. The limits on the i?-parity violation 



*The correction factor is less than 10% for bino masses under consideration [107 
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Figure 8: Upper hounds on the R-parity violation parameter derived from the Fermi LAT 
gamma-ray line limits in Fig. [3 For thermal leptogenesis, overproduction (O.P.) of gravitinos 
already excludes the left green region. As discussed in the text, the adopted strategy for line 
searches might fail for gravitino masses above a few hundred GeV (dashed black line). 



parameter ( are presented in Fig. El As discussed above, for gravitino masses above a few 
hundred GeV, our line search may fail due to a too large gamma-ray continuum contribu- 
tion, which is indicated by the dashed line. Note also that, at these high gravitino masses, 
the production of anti-protons in gauge boson fragmentation could further constrain the 
( parameter 35|, |36| . The implications of our gamma-ray line limits on the decay of the 



NLSP are discussed next. 



3.2.1 Stau NLSP 

In the case of the fi-NLSP, the total decay width of the lightest mass eigenstate is a mixture 
of left and right handed partial decays 



= sin^ 9f Tf + cos^ 9f T 



(19) 



Since the i?-parity breaking Yukawa couplings are typically proportional to the ordinary 
Yukawa couplings, decays into second and third familie s do minate. For definiteness, we will 
below assume a flavor structure as described in Ref. 107j . which is based on a Froggatt- 
Nielsen U(l) flavor symme try. In such a framework, the chiral state decays are dominated 
by the following channels 107 



The corresponding decay rates are 
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Figure 9: Contour plot of lower bounds on the stau NLSP decay length coming from the 
gamma-ray line constraints on the gravitino lifetime, as function of the stau and gravitino 
masses, m,fi and 7713/2 respectively. The lower gray region is excluded by electroweak precision 
tests (EWPT). For thermal leptogenesis, overproduction (O.P.) of gravitinos excludes at 
minimum the left green region, a limit which strengthens to the black-dashed line when 
assuming the universal boundary conditions (B), cf. Eq. ([T5]). The lower-right exclusion 
comes from the gravitino LSP requirement. Note that for high gravitino masses above a few 
hundred GeV our adopted line search strategy might overestimate the limits on the NLSP 
decay length, because the branching ratio into gamma-ray lines can become very small (see 
text). 



where the dimensionless parameter e is directly related to the i?-parity violating Yukawa 



couplings (see Ref. |107l | for details). In principle ( and e are independent parameters 



because they stem from different linear combinations of the bilinear i?-parity violating 
couplings fii, Bi and m^^. We choose ^ ~ e in order to set limits on the stau decay length. 
However, one has to keep in mind that in principle C, could be set to values much smaller 
than e by a proper choice of the parameters /ij, Bi and mf^. The behavior of the mixing 
angle Of with the stau mass can be deduced from the RGEs using the boundary 
conditions Eq. ( IT5l) . 

Using the upper limits on the i?-parity breaking parameter from Fig. [HI we can derive 
lower bounds on the stau decay length. Our results are shown in Fig. |9l The parameter 
space is already constrained by EWPT and overproduction bounds, and the lower limits 
on the neutralino decay length vary between 100 m and 10 km. It is interesting that if such 
particle were to be produced at the LHC, a sizable amount of their decay could take place in 



the detector [37[. We obtain a lowest possible decay length cTf ~ 85 m for 7713/2 — 16 GeV 



and mf^ ~ 100 GeV. 
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Figure 10: Like Fig. [PI but for a ncutralino NLSP. 



3.2.2 Neutralino NLSP 



A neutralino NLSP heavier than 100 GeV dominantly decays into W^i"^ and Z^i/ llOUlll . 
The corresponding decay width is directly proportional to the i?-parity breaking parameter 
C squared, which also enters the gravitino decay width in Eq. (ITH]) . As a consequence, the 



two quantities can be related through |107 



2^2 Ml G^M^m^o 2f{m^o,mw) + fim^o,mz)' 



where the phase space factor / is defined by 



,2 



/K,m,)={l-5 1 + 23) . (23) 



Using the gaugino mass relation from Eq. ( 1161) . lower bounds on the neutralino decay 
length cr^o can then be derived from Tab. [3] and Fig. |H Our results are summarized 
in Fig. [TOl For the parameter space allowed by EWPT and overproduction bounds, we 
obtain minimal decay lengths (9(100 m — 10 km), which are in the range of detectability 



of the LHC [37[. Decay lengths as small as cr-^o ~ 50 m are allowed for m^i2 — 17 GeV et 



m^o ~ 110 GeV. 



4 Conclusions 

We have presented a dedicated search for gamma-ray lines with energies between 1 and 
300 GeV in the diffuse gamma-ray fluxes measured by Fermi LAT between 4 Aug 2008 and 
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17 Nov 2010, based on the public 'DataClean' event class. Such lines could be produced in 
dark matter annihilation or decay, and they are highly relevant for dark matter searches 
since they provide an essentially background free signature. To calculate the significance 
of potential gamma-ray lines in the diffuse flux, we used a binned profile likelihood method 
and fitted the observed energy spectrum locally with a line signal as expected from decaying 
or annihilating dark matter inside the Galactic halo plus a power law that resembles the 
astrophysical background fluxes. No gamma-ray line was found at 5cr confidence level, and 
we derived conservative upper limits on the gamma-ray line flux that comes from regions 
with high Galactic latitudes or from the Galactic center (Tab. |3]). We also discussed the 
apparent presence of several line-like features at the few percent level in the public data, 
which appear at energies below 20 GeV with ~ 3a significance and which are most probably 
due to instrumental effects. 

Our conservative upper limits on the partial annihilation cross section of WIMP dark 
matter into 7-pairs, {av)^^^^^, improve the EGRET limits in the range = 1-10 GeV by 
up to an order of magnitude, while the previous Fermi LAT results {E^ = 30-200 GeV) are 
updated with two years of data and extended to the larger energy range E-y = 1-300 GeV 
(Fig. IS]). Despite the increased statistics and the choice of a sky-region with optimized 
signal-to-noise ratio, the previous Fermi LAT limits are however only slightly improved. 
The reason is a somewhat worse energy resolution of the public data with respect to the 
dedicated Fermi LAT analysis as well as our strategy to always quote conservative limits 
in certain energy bands. The lower bounds on the inverse decay widths of dark matter into 
monochromatic photons, T^_^^^, are at the level of 6 x 10^^ s over most of the considered 
energy regime (Fig. H]). Our results also constrain annihilation and decay channels with 
'-fZ^ final states. 

In a -R-parity breaking supersymmetric framework with gravitino LSP, we used the 
results of our gamma-ray line search to derive lower bounds on the gravitino lifetime. 
Considering supergravity models with two different types of universal boundary conditions 
at the grand unification scale, these lifetime constraints were used to set lower limits on 
the corresponding NLSP decay lengths. For gravitino and NLSP masses compatible with 
electroweak precision tests and overproduction constraints, we obtain (9(100 m— 10 km) 
lower bounds on the NLSPs decay lengths. Interestingly, all these decay lengths would be 
accessible at the LHC. 



Note. During the final stages of this work, first LHC limits on the CMSSM parameter 



space were published by the CMS collaboration [112[, which will further constrain the 
superparticle masses discussed in this paper. Furthermore, the anonymous referee pointed 
us to a talk about the ongoing gamma-ray line analysis by the Fermi LAT collaboration, 
where preliminary results are presented that are, where they overlap, in agreement with 
our findings, and where line-like artefacts at lower energies in the current public data are 
mentioned 



^See talk by E. Bloom at Aspen Winter Workshop on Indirect and Direct Detection of Dark Matter, 



Feb 2011, |http : //www. slac . Stanford. edu/exp/glast/aspenll/talks . asp 
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